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ABSTRACT 



We model the evolution of buoyant magnetic flux tubes in the Sun's con- 
vection zone. A flux tube is assumed to lie initially near the top of the stably 
stratified radiative core below the convection zone, but a segment of it is per- 
turbed into the convection zone by gradual heating and convective overshoot 

^ \ motions. The ends ( "footpoints" ) of the segment remain anchored at the base of 

r-| \ the convection zone, and if the segment is sufficiently long, it may be buoyantly 

unstable, rising through the convection zone in a short time. The length of the 

flux tube determines the ratio of buoyancy to magnetic tension: short loops of 

^ \ flux are arrested before reaching the top of the convection zone, while longer 

loops emerge to erupt through the photosphere. Using Spruit's convection zone 
model, we compute the minimum footpoint separation L^. required for erupting 
K^ \ flux tubes. We explore the dependence of L^ on the initial thermal state of the 

^r^ \ perturbed flux tube segment and on its initial magnetic field strength. Following 

an investigation of thermal diffusion time scales and the dynamic rise times of 

. _ , unstable fiux tube segments, we conclude that the most likely origin for magnetic 

O \ flux which erupts to the surface is from short length scale perturbations (L < L^) 

which are initially stable, but which are subsequently destabilized either by dif- 
fusion of heat into the tube or by stretching of the anchor points until L just 
exceeds Lc. In either case, the separation of the anchor points of the emergent 



H ' tube should lie between the critical distance for a tube in mechanical equilib- 

rium and one in thermal equilibrium. Finally, after comparing the dispersion of 
dynamic rise times with the much shorter observed active region formation time 
scales, we conclude that active regions form from the emergence of a single flux 
tube segment. 



1. Introduction 



I n several recent theories of solar and stellar dynamos (JGolub et al.lll98ll : iGalloway fc Weiss 



198ll : iDeluca fc Gilmarull986l ) magnetic flelds are generated not in the convection zone, but 
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in a thin, convectively stable layer just below it. In order for this magnetic flux to emerge to 
the photosphere where it can be observed, it must first make its way through the convection 
zone. It is generally believed that magnetic bu oyancy provides the force which pushes the 
magnetic flux toward the s urface ( Parkerl 1 19791 ) . Our goal in this paper is to use the thin 



flux tube approximation of ISpruitI ( 



198ll ) to study the emergence of buoyant loops of mag- 



netic flux whose ends are still anch ored in the stable layer below the convection zone. Our 
approach is similar to that used bv iMoreno-InsertisI (119861 ). Further d etails of our c alcula - 
tions may be found in IChou fc Fisheii (119891 ): henceforth Paper I, and iFisher et al.l ( 1l989l ). 
henceforth Paper II. Here we present only the essentials. 



2. Description of the Model 

In the thin flux tube approximation, the tube diameter is assumed to be smaller than 
any other relevant physical length scale, and the magnetic field is taken to be constant 
across the tube. Furthermore, since the magneto-acoustic transit time across the tube will 
be smaller than other time scales of interest, one may assume that there is always a balance 
between the gas pressure outside the tube and the sum of gas plus magnetic pressure inside 
the tube. One can think of a thin flux tube as a one- dimensional curve embedded in space. 
Our flux tube is tied down at both ends to the bottom of the convection zone, with the 
portion in between free to move in response to buoyancy, magnetic tension, and drag forces. 
We take the flux tube to lie in a vertical plane, defined by a horizontal x— axis at the base 
of the convection zone, and a vertical y— axis, and regard the height y of the the flux tube 
as a function of x and time t. The ends of the flux tube segment are anchored at x = and 
X = L. In all the calculations discussed in this paper, we assume that the flux tube initially 
{t = 0) lies fiat at the bottom of the convection zone {y = 0). In the discussion below, 
subscript e refers to the external plasma of the convection zone, and subscript i refers to the 
plasma inside the flux tube. 



2.1. The Basic Equations 

The velocity of the plasma in the flux tube consists of the speed of the tube itself in the 
direction normal to its length (denoted f„) and the speed of plasma moving along the flux 
tube (denoted by Vs)- The equation of motion for the velocity component Vn is given by 
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where the terms on the right hand side of equation ([T]) correspond, respectively, to magnetic 
tension, buoyancy, and aerodynamic drag forces. The drag coefficient C^ is taken to be 
unity. Note that the inertial term contains the "added mass" of th e displaced external 



plasma flowing around the flux tube (see, e.g. jLandau fc Lifshitall959l . §24). An insignicant 



centrifugal force due to parallel flows through bends in the flux tube has been neglected. 
Equation ([1]) must be supplemented by an equation relating the vertical motion of the tube 
to the velocity f„ and the flux tube slope dy/dx. From geometrical considerations, we flnd 
at a flxed value of x 

^ = v^^l + {dy/dxY . (2) 

The radius of curvature i?, the magnetic fleld strength _B, and the flux tube cross-sectional 
area A appearing in equation ([1]) are determined by 

1 d'^vldx'^ B^ 

_ _ <J y/ijx p _ p, _ i^. pA - ^ (Q\ 

R [l + idy/dxyf/'^' ' ' Stt' ' ^^ 

where Pi and Pe are the internal and external gas pressures, respectively, and $ is the total 
magnetic flux. The numerical techniques used for solving equations ([T]) and ([2]) are described 
in Paper I. 

The parallel flows of the plasma inside the flux tube can be calculated from equations 
(II]) and (l2l) and mass conservation arguments (see Paper I). These flows turn out to be 
sufficiently slow that the plasma can be taken to be in hydrostatic equilibrium. If we make 
the same assumption regarding the external plasma (i.e., ignoring effects of convection) , then 
it is straightforward to determine all the thermodynamic and magnetic variables as functions 
of height alone. This procedure is outlined below. It is not at all obvious a priori that the 
neglect of convective motions is justified. The presence of these motions could affect flux 
tube evolution in two ways. First, the convective motions could severely distort the shape 
of the flux tube. Second, convective motions introduce a turbulent pressure which perturbs 
hydrostatic equilibrium and this could affect the calculated magnetic fleld strength. We 
consider both effects in Paper I and conclude that the flrst effect could be important for 
some of the thinnest flux tubes we study, whereas the second is not important for any of the 
cases we have studied. 

The structu r e of th e background external atmosphere is taken to be the convection zone 



model of ISpruitI (J1974J ). This atmosphere is slightly superadiabatic due to convective heat 
transport. The information in Spruit's paper is sufficient to calculate the thermodynamic 
and magnetic quantities needed for our dynamic model. We assume that the plasma inside 
the flux tube is isolated from heat transfer and behaves adiabatically. There will of course 
be heat transport between the flux tube and its surroundings, but we have not included 
this specifically in the model. Our approach in this paper is to approximate the effects of 
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heat conduction by choosing a physically self-consistent value of the temperature difference 
between the flux tube and its environs {i.e. the parameter 77 introduced below) at the base of 
the convection zone. This will be described further in the section on time scales for emerging 
flux tubes. 

It is straightforward to derive an equation for /3 [/3 = SiiPe/B'^ = Pe/{{Pe — Pi)] in 
terms of the height variation of the external temperature difference 5T = Tg — Ti between 
the external and internal plasma; this avoids the numerical difficulty of subtracting two gas 
pressures which are nearly equal. Starting from the equilibrium equations 

-f- = -Pe9\ -1- = -Pi9 4 

dy dy 

one finds 

1 df3 1 ST/Te 



(5) 



I3{I3 -l)dy A 1 - 5T/T, 
where A is the external pressure scale height. This has the solution 

m'' = /3o"'exp[-7(y)] + {1 - exp[-7(y)]} (6) 

where /3o = /3(0) and 

™' h A(!/0 1 - my')IUy') ■ * ' 



The temperature difference is simply 



y 



6T{y) = 6To- dy\V - Vad)f^9/Rg (8) 

Jo 

where 6T0 = ST{0), V = (dlriT/dlnP), Vad = {dlriT / dlnP) s , and /i and Rg are the mean 

mass per particle and the gas constant. The temperature difference at the bottom of the 

convection zone STq is given in terms of parameters rj and Bq by 

6To = vTe{0)/f3o, (9) 

where /3o = SnPe/BQ. The quantity Bq is the magnetic field strength at the base of the 
convection zone, and is one of the free parameters of the problem, while 77 is a dimensionless 
measure of the temperature difference at the base of the convection zone, and is constrained 
to lie between and 1. When 77 = 0, the internal temperature is equal to that of the 
surrounding plasma, whereas for rj = 1, the density inside the fiux tube is equal to that 
outside. Together with L and $, specification of Bq and rj completely defines the fiux 
emergence problem in our model. 

Values of all these parameters are highly uncertain, so we have attempted to cover a 
wide range of parameter space. For most of our simulations, Bq was chosen to be 10*^ G, 
10^ G, or 10^ G, and the fiux $ was taken as 10^^ Mx, lO^o Mx, or 10^2 Mx. Our choice of 
these values is discussed further in Paper I. 
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2.2. Example of a Simulation 

To demonstrate our model, we briefly discuss the results of a single simulation. In this 
case, Bo = 10^ G, $ = 10^® Mx, r/ = 0, and L = 2 x 10^° cm. The position of the tube at 
numerous times is shown in Figure [1] This flux tube is "unstable" and emerges from the top 
of the calculational domain without reaching an equilibrium configuration. Note that the 
upward motion of the flux tube slows significantly when its apex is roughly half-way through 
the convection zone. This indicates that the magnetic tension was almost strong enough to 
balance the buoyancy force. Had L in fact been only slightly smaller (< 1.97 x 10^^ cm), 
the flux tube would have been "stable" and would have stopped rising when the height of 
the flux tube apex y^ ^ 9 x 10^ cm. We explore this interesting behavioral dichotomy in the 
following subsection. 



3. The Critical Length Scale 

From the discussion in the previous section, it is apparent that for given values of the 
other parameters, a critical length L^ exists. For L < Lc, the flux tube reaches a stable 
equilibrium, while for L > L^, an unstable eruption to the photosphere occurs. Therefore, a 
knowledge of the dependence of Lc on the other parameters is very important for understand- 
ing the nature of flux emergence. The magnetic flux $ affects the drag force and therefore 
the rate of rise, but it does not affect the relative balance between buoyancy and tension, so 
Lc is independent of $. To determine the dependence of Lc on rj and Bq, we have used the 
numerical model described in the previous section. In Figure [21 we plot Lc as a function of 
T] for Bo = 10^ G and Bq = 10^ G. 



c 



There are several features of Figure [2] worthy of mention. First, note that for fixed Bo, L 
increases as rj increases from to 1. This is easily explained by the decrease in the relative 
strength of magnetic buoyancy compared to magnetic tension as the temperature inside 
the flux tube decreases. Thus the footpoint separation L over which magnetic tension can 
overcome buoyancy is increased (recall that for fixed apex height tension will scale roughly 
as I/"^). Second, as t] approaches 1 for fixed Bo, Lc approaches a finite value, even though 
the buoyant force goes to zero at the base of the convection zone. Since for r] = 1, the initial 
horizontal state is an equlibri um configuration, one can perfo rm a linear stability analysis 



of equations ([T]) and ([2]) (see ISpruit fc van Ballegooijenl Il982l . and Paper II) to obtain Lc 
analytically 

Lc = v/2/(l-V)7rA , (10) 

which is plotted as the horizontal long-dashed line in Figure [2J This analytic result corre- 
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Fig. 1. — Plot of the flux tube shape at times separated by 2 x 10^ s. 
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sponds very well to our numerical result for r] = 1 and Bq = 10 G, but it is about 10% 







smaller than the critical length (for emergence of the flux tube to the photosphere) which 
we find numerically for rj = 1 and Bq = 10^ G. The reason for this interesting discrepancy 
is discussed in Paper II. 



4. Time Scales for Emerging Flux Tubes 

If magnetic flux emerges from the convective overshoot region and evolves according 
to our model, only those flux tubes we have labeled "unstable" [L > L^) can reach the 
photosphere. The time required for the flux tube to rise through the convection zone is 
therefore of great interest. We assume that a flux tube is formed in the stably stratified region 
beneath the convection zone, and begins its buoyant rise when perturbed upward into the 
convection zone by gradual heating and convective overshoot motions. Before perturbation, 
the flux tube must be in mech anical equilib rium (pe = Pi) as otherwise it would adjust by 



rising or sinking adiabatically ( lParkeiill979l ). The nature of a perturbation is described by 
the two parameters L and r]. L is the length scale of the perturbation, and t] measures its 
thermodynamic state. When rj = 1, the perturbed flux tube is still in mechanical equilibrium 
{pi = Pe)i whereas if 77 = 0, the flux tube is in thermal equilibrium with its surroundings 
(Tj = Te). The proper choice of r] for the perturbation can be made by comparing the 
thermal diffusion time at the bottom of the convection zone Tth to the dynamic rise time 
Td for unstable 77 = 1 flux tubes. If Tth >> ^d, then r] = 1 is the best description of the 
perturbation, whereas if Tth « Td, then r^ = is most appropriate. When the two are 
comparable, an intermediate value of r] is called for. The thermal diffusion time scale is 
given by (see Paper II) 

^ ^^^Picy^n 

ti^ 16 aT^ ' ^ ' 

where a is the Stefan-Boltzmann constant, kr is the Rosseland mean opacity [cm^ g~^], and 
a = [^/{71'Bq)Y^'^ is the radius of the flux tube at the base of the convection zone. The 
dynamic rise time r^ for r] = 1 can be obtained from the numerical simulations. Note that 
the thinner (smaller $) flux tubes will have shorter thermal time scales and longer dynamic 
time scales. We find that Td does not depend strongly on L, provided L exceeds Lc by at 
least 5%. We therefore have used L = 4 x 10^° cm (cf. Figure |2]) to compute our estimates 
of Td- Computed values of r^ and Tth as functions of $ and Bq are shown in Table 1, along 
with the resultant values of rj. We also include in Table 1 the "corrected" dynamical rise 
times (denoted r^) using these self-consistent values of rj in the numerical simulations. 

From the information in Table 1 and the critical length scales of Figure El it is now 
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Fig. 2. — Critical length scale as a function of rj for Bq = 10^ and 10^ G. The long- 
dashed horizontal line corresponds to the value from equation (TTOl) . and the short-dashed 
line corresponds to Parker's (1979) isothermal result of 27rA. 
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possible to construct a general picture of flux tube evolution. For given values of $ and Bq, 
perturbation length scales divide naturally into three ranges: short, L < Lc{r] = 0); medium, 
Lciv = 0) < L < Lc{rj = 1); and long, Lc{ri = 1) < L. If Tfh « t^ then only short flux 
tubes are stable, while longer flux tubes rise to the photosphere on the dynamical time scale. 
But if Tth » Td, both short and medium flux tubes are stable. Only the short flux tubes are 
"absolutely" stable, however, as the medium length tubes are able to rise quasi-statically on 
the thermal time scale as heat leaks into them. This is equivalent to gradually decreasing t] 
from 1 to 0. However, before t] reaches zero, a medium-L flux tube will suddenly find itself 
with L > Lc, and will then erupt to the photosphere on the much shorter dynamic rise time. 

In the most straightforward {e.g. au or a'^u) kinematic dynamo models the oscillatory 
nature of the solar cycle is due to periodic conversion of poloidal to toroidal magnetic field 
and back via the mechanisms of stretching by differential rotation and the a-effect (which 
is usually attributed to net helicity in convective or convective overshoot motions). The 
dynamo period is therefore at least as great as the time necessary to reorient and stretch 
poloidal field into toroidal field. For these dynamos to work, therefore, magnetic flux must 
remain in the dynamo region long enough for reorientation to occur, i.e., for a significant 
fraction of the solar cycle. The conclusions drawn below are based on this premise. An 
alternative picture is that the dynamo generates fields on a much shorter time scale, in 
which case our arguments below do not hold. In that case, however, the period of the solar 
cycle itself remains unexplained and some additional unknown mechanism must be invoked 
to account for it. 

From a comparison of the time scales in Table 1 with the length of half a solar cycle 
(roughly 11 years or 3.4 x 10^ s), some tentative conclusions about the nature of the flux 
tube perturbations can be drawn. First, note that the rise time r^ of unstable flux tubes 
is always much shorter than the duration of the solar cycle. If there were a continuous 
source of perturbations with length scales L > L^, then magnetic flux could not remain 
stably submerged on time scales much longer than r^. One therefore concludes that such 
perturbations are either rare or nonexistent. Indeed, if convective overshoot motions are 
responsible for the flux tube perturbations, one would expect typical perturbation length 
scales to be of order the eddy size, which in mixing length theory is roughly the pressure 
scale height near the bottom of the convection zone, i.e. 6 x 10^ cm. This is significantly less 
than any of the Lc values in Figure |2l A perturbation with L > Lc would require coherent 
action on the part of several adjacent eddies. Assuming this to be unlikely, the magnetic flux 
which does eventually emerge must originate from perturbations in the short or intermediate 
L range. We consider two possibilities. 

First, if the separation L between anchor points is truly fixed in time, then only those 
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perturbations in the intermediate L range can emerge at all, since short L perturbations 
result in completely stable structures. Furthermore, only those combinations of Bq and $ 
which give rise to 77 = 1 in Table 1 are viable, i.e. those flux tubes with $ k. 10^^ Mx. In 
that case, we expect magnetic flux to emerge to the surface on thermal diffusion time scales 
and to have a footpoint separation between Lc{r] = 0) and Lc{ri = 1). 

The second possibility is that there is stretching of the anchor point separation L past 
Lc after a stable short L perturbation has been made. There are two possible mechanisms 
for this. In the first instance, stretching could be accomplished by differential rotation, 
for example, if the anchor points differ in latitude. The second mechanism is the thermal 
heating of the anchor points themselves. This may cause a gradual "unzipping" of the 
flux tube segment, until the footpoint separation exceeds L^. The time scale for this to 
occur depends on a number of factors, including how deeply buried the remainder of the 
flux tube is in the overshoot region. This is discussed further in Paper II. Here again, for 
either mechanism, we expect the anchor point separation of emerging flux to lie in the range 
L,{7] = 0) < L < L,{r] = 1). 

As a final point, note that the dispersion in rise times for individual flux tubes with 
different values of $ (a month or more - see Table 1) is much greater than the formation 
times of active regions (typically 2-3 days). We believe this indicates that active regions are 
formed from a single emerging flux tube. If many tubes of differing sizes were perturbed 
simultaneously, or if a single flux tube became highly fragmented near the bottom of the 
convection zone, one would expect the emergence of flux to be dispersed over a much greater 
time corresponding to the difference in rise times between small $ {e.g., 10^^ Mx) and large $ 
{e.g., 10^^ Mx) flux tubes. Any fragmentation therefore probably takes place in the topmost 
portion of the convection zone when the tube is emerging most rapidly. 



5. Conclusions 



We have developed a model, based on the thin flux tube approximation of lSpruitI ( 119811 ). 
for studying the emergence of magnetic flux through the convection zone when the footpoints 
of the flux loop are anchored a distance L apart in the stable layers below. Figure [1] shows 
the evolution of a flux tube whose footpoints are sufficiently far apart {L > L^) to allow it 
to rise to the top of the convection zone. We have used our model to explore the critical 
length Lc separating flux tubes which form stable magnetic loops in the convection zone 
{L < Lc) from those which erupt through the photosphere {L > Lc). Figure [2] shows the 
dependence of Lc on the temperature defect parameter rj for two values of the field strength 
at the base of the convection zone Bq, 10^ G, and 10^ G. We find that for perturbations on 
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Table 1: Flux Tube Time Scales 





5o = 10^ [G] 


Bq = 10^ [G] 


Bo = W [G] 


Uncorrected (77 = 1) rise t 


ime T(i [s 






$ = 10^2 [Mx] 


8.3 X 10^ 


6.4 X 10^ 


3.9 X 10^ 


$ = 10^0 [Mx 


1.8 X 10^ 


1.1 X 10^ 


4.7 X 10^ 


$ = lO^s [Mx 


5.3 X 10^ 


3.0 X 10^ 


8.1 X 10^ 


Thermal diffusion time Tty 


I [^. 


$ = 1022 [Mx 


2.7 X 10^ 


2.7 X 10^ 


2.7 X 10^ 


$ = 1020 [Mx 


2.7 X 10^ 


2.7 X 10^ 


2.7 X 10^ 


$ = lO^s [Mx 


2.7 X 10^ 


2.7 X 10^ 


2.7 X 10^ 


Self-consistent t] 


$ = 1022 [Mx 


1.0 


1.0 


1.0 


$ = 1020 [Mx 


0.0 


0.0 


0.5 


$ = lO^s [Mx 


0.0 


0.0 


0.0 


Corrected rise time t^. [s 




$ = 1022 [Mx 


8.3 X 10^ 


6.4 X 10^ 


3.9 X 10^ 


$ = 1020 [Mx 


6.7 X 10^ 


3.1 X 10^ 


8.6 X 10^ 


$ = lO^s [Mx 


2.1 X 10^ 


9.8 X 10^ 


1.7 X 10^ 
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length scales L > Lc{Bo,ri) the rise time is short compared to the duration of the solar cycle. 
Under the assumption that the solar cycle can be modeled as a kinematic {e.g., au or a'^u) 
dynamo model operating below the convection zone, we conclude that such perturbations 
must therefore be rare or nonexistent. We speculate that the magnetic flux which does erupt 
through the photosphere forms initially from perturbations with L < Lc, resulting in stable 
structures. These are subsequently destabilized either by thermal diffusion or by stretching 
of the anchor points until L exceeds L^. In either case, we expect that the anchor point 
separation L should fall roughly within the range Lc{i] = 0) to Lc{i] = 1). The rise time of 
flux tubes is the time scale for conduction of heat into the tube in the flrst case, and the 
stretching time in the second case. Finally, we argue that active regions are formed from the 
emergence of a single flux tube segment. 
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